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Characteristics  of  lean  premixed  conventional  micro-combustion  and  lean  non-premixed  flameless 
regime  of  methane/air  are  investigated  in  this  paper  by  solving  three-dimensional  governing  equations. 
At  moderate  equivalence  ratio  (0  =  0.5),  standard  k-E  and  the  eddy-dissipation  concept  are  employed  to 
simulate  temperature  distribution  and  combustion  stability  of  these  models.  The  effect  of  bluff-body  on 
the  temperature  distribution  of  both  conventional  and  flameless  mode  is  developed.  The  results  show 
that  in  the  premixed  conventional  micro-combustion  the  stability  of  the  flame  is  increased  when  a  trian¬ 
gular  bluff-body  is  applied.  Moreover,  micro-flameless  combustion  is  more  stable  when  bluff-body  is 
used.  Micro-flameless  mode  with  bluff-body  and  7%  02  concentration  (when  N2  is  used  as  diluent)  illus¬ 
trated  better  performance  than  other  cases.  The  maximum  temperature  in  premixed  conventional  micro¬ 
combustion  and  micro-flameless  combustion  was  recorded  2200  K  and  1520  K  respectively.  Indeed,  the 
flue  gas  temperature  of  conventional  mode  and  flameless  combustion  was  1300  I<  and  1500  K  respec¬ 
tively.  The  fluctuation  of  temperature  in  the  conventional  micro-combustor  wall  has  negative  effects 
on  the  combustor  and  reduces  the  lifetime  of  micro-combustor.  However,  in  the  micro-flameless  mode, 
the  wall  temperature  is  moderate  and  uniform.  The  rate  of  fuel-oxidizer  consumption  in  micro-flameless 
mode  takes  longer  time  and  the  period  of  cylinders  recharging  is  prolonged. 
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1.  Introduction 

3.1.  Development  of  micro-scale  combustion 

Recently,  the  demand  for  lightweight  and  compact  energy 
resources  has  significantly  increased  because  high  energy  density 
could  not  be  supplied  by  traditional  batteries.  The  demand  for  com¬ 
pact  energy  is  expected  to  be  increased  due  to  development  of  elec¬ 
tronic  devices  and  thus  micro-power  will  be  required  to  supply 
energy  of  these  enhanced  functionalities  [1,2].  Since  the  energy  den¬ 
sity  (per  unit  mass  or  unit  volume)  of  hydrocarbon  fuels  is  much 
greater  than  the  traditional  batteries,  micro-power  generation 
based  on  combustion  of  hydrocarbons  has  been  more  considered. 
Furthermore,  the  negative  impacts  of  conventional  batteries  on 
the  environment  upon  disposal  highlight  the  necessity  of  emergence 
of  new  micro-power  generation  methods  for  micro-electro- 
mechanic-systems  (MEMS)  [3],  Although,  the  package  of  combus¬ 
tion  based  micro-power  generation  could  be  significantly  shrank, 
the  demanded  power  is  not  compromised  if  the  chemical  energy 
of  the  fuel  is  utilized  efficiently.  Despite  the  low  efficiency  of 
micro-combustion  systems,  the  priorities  of  combustion  based 
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micro-power  generating  in  a  few  Watts  within  an  extremely  small 
volume  was  noticed  [4],  In  order  to  address  the  raising  demand  for 
small-scale  electricity  sources,  various  micro-combustion-based 
power  generators  have  been  developed.  Hence,  enhancement  of 
flame  stability  and  thermal  efficiency  of  micro-combustors  has 
become  new  challenge  in  combustion  investigations  [5,6],  The  fun¬ 
damental  concepts  of  combustion  characteristics  of  micro-scale 
combustors  play  crucial  role  in  the  improvement  of  system  effi¬ 
ciency  and  optimization  of  the  design.  Since  combustor  volume  is 
reduced  in  micro-scale  combustion,  conspicuous  radical  destruction 
and  heat  loss  from  micro-combustor  walls  are  expected. 
Accordingly,  application  of  micro-thermophotovoltaic  (TPV)  power 
generation  in  MEMS  has  been  enhanced.  Elimination  of  moving 
parts,  high  reliability  and  highly  robust  of  TPV  are  the  main  charac¬ 
teristics  of  TPV  which  makes  it  suitable  for  application  in  commer¬ 
cial  electronic  devices  [7].  In  the  recent  years,  conspicuous 
progress  has  been  obtained  in  micro-scale  combustion  experimen¬ 
tally  and  numerically.  These  investigations  have  helped  to  under¬ 
stand  the  fundamentals  of  micro-combustion  in  terms  of 
flammability  limits,  flame  stability,  emitter  and  thermal  efficiency. 
Application  of  micro-combustion  systems  in  MEMS  was  successfully 
experimented  by  Waitz  et  al.  [8]  in  micro-turbines  and  Yang  et  al.  [9] 
in  the  TPV  in  the  last  decade.  Since  the  ratio  of  surface  to  the  volume 
of  micro-combustor  is  higher  than  regular  combustors,  the 
possibility  of  thermal  quenching  in  micro-combustors  is  very  high 
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Nomenclature 

Pt 

turbulent  viscosity 

ct 

time  scale  constant  (=  0.4082) 

Pk 

turbulence  kinetic  energy  production 

turbulent  Prandtl  numbers  for  k 

Pb 

effect  of  buoyancy 

Oe 

turbulent  Prandtl  numbers 

Ym 

the  contribution  of  the  fluctuating  dilatation  to  the 

Sfc.  se 

user-defined  source  terms 

overall  dissipation  rate 

Cie.  c2e,  C3e  constants 

Ts,o 

outer  surface  temperature 

Gk 

the  generation  of  turbulence  kinetic  energy  due  to  the 

Too 

ambient  temperature  set  at  300  K 

mean  velocity  gradients 

h 

natural  convection  coefficient  considered  constant  va- 

Gb 

the  generation  of  turbulence  kinetic  energy  due  to 

lue  5  W/m2  K 

buoyancy 

G 

Stefan-Boltzmann  constant  (=  5.67  ■  10~8  W/m2  K4) 

7] 

diffusion  flux  of  species  j 

e 

solid  surface  emissivity 

l<eff 

effective  conductivity 

Q 

volume  fraction  constant  (=  2.1377) 

Ym 

the  contribution  of  the  fluctuating  dilatation  in  com- 

V 

kinematic  viscosity 

pressible  turbulence  to  the  overall  dissipation  rate 

due  to  the  extremely  high  heat  loss  from  flame  [10].  Therefore,  sta¬ 
bility  of  the  flame  in  combustors  with  sub-millimeter  size  has 
become  a  new  challenge  [11],  Indeed,  the  stability  of  the  flame  in 
micro-scale  combustors  is  influenced  by  absorption  and  destruction 
of  generated  combustion  radicals  [12  .  Hence,  feasibility  of  using 
catalyzed  combustion  [13]  as  well  as  application  of  external  heat 
[14  and  heat  recirculation  [15]  were  experimented.  It  was  found 
that  flame  stability  in  micro-scale  combustors  depends  strangely 
on  the  heat  recirculation  through  combustor  walls  [16],  Preheating 
the  reactants  via  heat  transfer  from  high  temperature  region  broads 
the  reaction  zone  in  micro-scale  combustion  [17”.  To  decrease  heat 
loss  to  the  ambient  and  impalement  heat  recirculation  to  the 
unburned  zone,  specific  material  should  be  selected.  Therefore, 
material  selection  in  micro-scale  combustors  is  vital  due  to  strong 
thermal  coupling  between  reacting  mixture  and  micro-combustor 
walls.  Maruta  et  al.  [18]  pointed  out  that  flame  instabilities  in  non- 
premixed  micro-scale  combustion  could  be  attributed  to  the  inter¬ 
action  between  flame  structure  and  flame  flow,  heat  loss  and  mass 
transfer  limitations.  Flame  stability  and  characteristics  of  non-pre- 
mixed  sub-millimeter  combustor  was  investigated  by  Prakash 
et  al.  [19],  The  source  of  flame  instabilities  in  premixed  combustion 
of  narrow  channels  was  investigated  by  Sanchez-Sanz  et  al.  [20], 
Since  the  size  of  combustor  shrinks  to  sub-millimeter  in  micro¬ 
scales,  some  combustion  parameters  such  as  flame  ability,  flame  sta¬ 
bility,  flame  thickness,  heat  loss,  flow  field  and  thermal  efficiency  are 
changed.  Therefore,  the  role  of  numerical  and  computational  inves¬ 
tigation  has  been  more  highlighted  in  micro-combustion  due  to  dif¬ 
ficulty  of  measuring  various  quantities  in  such  limited  circumstance 
in  experimental  investigations.  The  impacts  of  various  parameters 
such  as  combustor  size,  the  geometry  of  combustor  and  boundary 
conditions  on  the  flame  temperature  of  hydrogen-air  premixed 
micro-combustion  were  investigated  by  Li  et  al.  [21]  numerically. 
The  authors  stipulated  that  flame  temperature  in  micro-scale  com¬ 
bustion  could  be  influenced  by  the  mixture  flow  rate,  the  size  and 
the  geometry  configuration  of  micro-combustor.  Karagiannidis 
et  al.  [22]  confirmed  that  a  wide  range  of  flammability  could  be 
offered  in  catalytic  micro-combustors.  The  effects  of  equivalence 
ratio  as  well  as  inlet  velocity  on  the  structure  of  the  premixed  CH4/ 
air  flame  in  micro-scale  size  were  investigated  by  Feng  et  al.  [23] 
computationally.  It  was  found  that  the  position  of  the  flame  is 
shifted  to  the  downstream  in  high  inlet  velocities.  Indeed,  the  high¬ 
est  flame  temperature  could  be  achieved  in  stoichiometric 
condition. 

3.2.  Flameless  combustion  technology 

The  commencement  of  flameless  combustion  research  traced 
back  to  about  1990s,  mainly  attributed  to  the  augmentation  of 


the  environmental  concerns.  Depletion  of  fossil  fuel  resources 
and  anxieties  about  environmental  issues  conducted  combustion 
research  community  to  find  environmentally  friendly  combustion 
technology  and  optimized  fuel  consumption  [24],  Flameless  com¬ 
bustion  [25]  or  so-called  moderate  intensive  low  oxygen  dilution 
(MILD)  combustion  [26]  or  high  temperature  air  combustion  tech¬ 
nology  (HITAC)  [27]  has  been  developed  in  the  resent  years 
because  of  its  ability  to  combine  high  performance  of  combustion 
with  extremely  low  pollutant  formation  [28],  The  uniform  temper¬ 
ature  inside  the  chamber  emerges  perfectly  stirred  reactor  circum¬ 
stances.  Since  the  peak  of  temperature  reduces  in  flameless  mode 
and  hot  spots  are  eliminated,  NOx  formation  reduces  dramatically 
[29],  Indeed,  low  emission  of  flameless  mode  is  attributed  to  the  air 
dilution  by  inert  gases  such  as  C02  and  N2  [30,31],  In  flameless 
combustion,  the  temperature  of  inlet  diluted  oxidizer  is  higher 
than  self-ignition  of  the  fuel,  therefore  ignition  is  eliminated  [32], 
Due  to  the  significant  potentials  of  the  flameless  mode  to  deal  with 
various  fuels,  implementation  of  flameless  combustion  in  micro¬ 
scale  power  generation  could  solve  dilemma  problems  in  micro¬ 
combustions  and  thus  application  of  this  model  of  combustion  in 
micro-scale  combustors  could  enhance  combustion  stability  and 
thermal  efficiency  [33,34],  In  the  flameless  mode,  the  reaction  zone 
is  distributed  throughout  the  chamber  consequently;  the  fluctua¬ 
tions  of  temperature  are  eliminated  and  uniform  low  temperature 
is  observed.  Furthermore,  low  oxygen  concentration,  elimination  of 
audible  and  visible  flame,  low  level  of  soot  formation  are  the  main 
characteristics  of  flameless  mode  [35],  Sine  flameless  mode  shows 
complicated  performance,  analytical  and  computational  investiga¬ 
tion  of  this  regime  has  received  more  attention  especially  in 
macro-scale  combustors.  Kim  et  al.  [36]  compared  different  reac¬ 
tion  mechanisms  in  natural  gas  flameless  combustion  by  using  four 
different  global  reaction  mechanisms  and  setting  standard  k-e 
model  and  the  eddy-dissipation  concept  (EDC).  Dally  [37]  stipu¬ 
lated  that  based  on  the  numerical  results  of  centerline  tempera¬ 
tures  of  the  combustor  achieved  by  adoption  of  k-e  and  flamelet 
model,  this  settings  are  in  good  agreement  with  experimental 
results  for  CH4  flameless  combustion.  Christo  and  Dally  et  al.  [38] 
pointed  out  that  it  is  necessary  to  consider  differential  diffusion 
effects  in  CFD  modeling  due  to  of  calculations  accuracy  and  the 
numerical  results  of  eddy  dissipation  concept  (EDC)  solver  shows 
higher  accuracy.  Hosseini  and  Wahid  [39]  investigated  various 
aspects  of  biogas  flameless  combustion  numerically.  It  was  found 
that  reaction  time  and  Damkohler  number  plays  crucial  role  in 
the  computational  settings.  Although,  characteristics  of  flameless 
combustion  in  macro-scale  combustors  have  been  developed 
numerically  and  experimentally,  few  documents  could  be  found 
about  micro-scale  flameless  mode.  In  the  present  paper, 
micro-flameless  formation,  temperature  distribution  in  the 
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micro-flameless  combustor,  the  effects  of  bluff-body  on  the  perfor¬ 
mance  of  micro-flameless  mode  and  the  effects  of  diluent  on  the 
micro-flameless  formation  are  investigated  numerically  and  the 
results  are  compared  to  premixed  micro-combustion. 

2.  Methodology 


treat  as  continuums  and  thus  Navier-Stokes  equations  are  valid 
[42].  Pervious  investigations  indicated  that  the  stability  of  turbu¬ 
lent  micro-combustion  is  higher  than  laminar  one  [43,44],  In  addi¬ 
tion,  enhancement  of  turbulence  aids  flameless  regime  to  be 
sustained.  Therefore,  standard  k-e  turbulence  model  with  follow¬ 
ing  transport  equations  is  applied  to  model  flameless  micro-com¬ 
bustion  system  [39], 


A  three-dimensional  (3D)  micro-combustion  is  simulated  by 
ANSYS  14  using  ANSYS  Modeler  to  design  the  micro-scale  chamber 
and  ANSYS  Meshing  to  mesh  the  combustor  [40].  In  order  to 
improve  the  convergence  rate  and  scalar  properties,  mesh  refine¬ 
ment  is  employed  and  therefore,  grid  resolution  for  smooth  flow 
representation  is  ensured.  In  this  investigation,  a  circular  tube  with 
2  mm  diameter  is  considered  as  a  non-premixed  micro-combustor. 
It  is  assumed  that  fuel  jet  is  surrounded  by  the  oxidizer,  the  diam¬ 
eter  of  inlet  fuel  is  1  mm  and  the  length  of  micro-combustor  is 
15  mm.  To  simulate  premixed  conventional  micro-combustor,  a 
tube  shape  combustor  with  2  mm  diameter  is  considered  where 
reactants  are  injected  from  the  inlet  step  with  1  mm  diameter 
because  it  was  demonstrated  that  performance  of  conventional 
micro-combustion  could  be  intensified  by  using  inlet  step  [41], 
Fig.  1  illustrates  the  schematic  of  studied  micro-combustors  and 
its  unstructured  mesh. 

In  this  steady-state  CFD  simulation,  methane  (CH4)  is  assumed 
as  the  fuel  and  the  swirl  velocity,  pressure  work  and  viscous  forces 
are  neglected.  The  equivalence  ratio  (0)  of  both  premixed  conven¬ 
tional  combustion  and  non-premixed  flameless  mode  is  0.5.  Since 
the  length  of  micro-combustor  is  longer  than  the  molecular  mean- 
free  path  of  the  reactants  through  the  micro-combustor,  reactants 


3,  d  ,  .  .  d 
s(pk)  +  fi*(pkU')  =  sr 


dt 


uc\  dk 

A <  + 


Ok)  dxj 


dx< 


-  pe  -  YM  +  Sk 


+  Pk  +  Pb 


(1) 


d  .  d  d 

^AP£)+^iP£U>)  = 


dt 


dXi 


dXj 


P  ~F 


llt\  Pk 


aj  dXj 


+  Ci  e^(Pk  +  C3ePb) 


~  C2ep-r+  Se 


(2) 


Sk,  Se  user-defined  source  terms,  ClEI  C2e,  C3e  are  constants.  Where 
turbulent  viscosity  and  turbulence  kinetic  energy  production  are 
modeled  by  Eqs.  (3)  and  (4)  respectively. 


r  1(2 

Pt=PCMJ 


(3) 


Pk  =  -pu'u’ 


dUj 

dUj 


(4) 


The  accuracy  and  robustness  of  k-e  model  for  turbulent  flows  is 
very  high.  Thus,  the  accuracy  of  CFD  modeling  could  be  enhanced 
due  to  fully  turbulent  circumstance  in  flameless  combustion 
chamber  [45], 
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3D  wire  frame  figure  of  the  inlet 
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Fig.  1.  The  schematic  of  micro-combustor  (A)  non-premixed  combustor,  (B)  premixed  combustor  and  (C)  unstructured  mesh  of  the  micro-combustor. 
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Eqs.  (5)  and  (6)  named  Reynolds-average  Navier-Stokes  (RANS) 
equations  are  employed  for  continuity  and  momentum  in  turbu¬ 
lent  condition. 


Reaction  rate  is  calculated  by  Arrhenius  equation  [47], 
k  =  AT''e 


(15) 


dp 

dt 


+ 


d_ 

dXi 


(put)  =  0 


(5) 


d_ 

dt 


(PU<) 


d_ 

dXj 


(pu.uj) 


dp  d  I-  f&Ui  dUj  2  dUj\ 
Wi  +  dxj[fi\dxj  +  dxi~3diidXi) 


(6) 


The  energy  equation  in  turbulent  circumstance  developed  by 
Ref  [39]  is  applied  in  this  case. 

|  (pE)  +  V  •  (V(PE  +  p))  =  V  ■  (keffVT  -  Jjtfi  +  (Teff  ■  V))  +  Sh 


(7) 

where  keff  is  effective  conductivity,^  is  diffusion  flux  of  species  j,  Sh 
is  a  combination  of  heat  chemical  reactant  and  every  other  heat 
resource  which  are  defined. 

On  the  right  hand  side  of  Eq.  (7),  the  first  three  terms  represent 
energy  transfer  due  to  conduction. 

Second-order  scheme  is  set  to  solve  related  governing  equa¬ 
tions.  The  residual  of  all  variables  is  adopted  1CT3  except  energy 
which  is  set  10~6  to  ensure  the  convergence  of  the  solution.  The 
swirl  velocity  is  neglected  in  both  conventional  flame  and  flame¬ 
less  mode.  Since  no  remarkable  difference  in  terms  of  flame  tem¬ 
perature  was  reported  in  fully  developed  and  uniform  inlet 
velocity  of  the  mixture  in  the  micro-combustor  [21],  the  inlet 
velocity  of  mixture  is  considered  uniform. 

Methane-air-2step  by  following  equations  [46]  is  employed  to 
model  species  transport. 

CH4  +  3/2O2  -  C0  +  2H20  AG  = -632.68  kj/mol  (8) 

CO  +  I/2O2  C02  AG  =  -257.19  kj/mol  (9) 

The  variation  of  equilibrium  constant  against  the  temperature 
can  be  written  as  Eq.  (10). 


I<  = 


rxw 

iLwr 


(10) 


where  Xt,  XJt  vt  and  Vj  are  the  molar  fractions  of  the  reactant  i,  prod¬ 
uct  j,  the  stoichiometric  coefficients  of  the  reactant  i  and  product  j 
respectively. 

For  computation  of  Arrhenius  parameters  of  reaction  (9)  (as  an 
example)  the  equilibrium  constant  Ki  and  I<2  are  defined  as: 


(XCo2)^  Ko2 

(Xo2)^(Xco)I'“  X°25X’0 


(11) 


I<2 


XU.D  yl 

02  Aco 


X 


1 

co2 


1 

X, 


(12) 


The  backward  rates  and  the  forward  rates  are  derived  from 
equilibrium  constants  of  reaction  as  shown  in  Eq.  (13). 


k 


(13) 


( b  for  backward  and  /  for  forward  reaction). 
Therefore, 


kft  (^[02]1/2 [CO]1  -^-[C02]'j  =  kf2k,  ^[02]1/2[CO]1  -^-[C02]1^ 

(14) 


where  k  is  reaction  rate,  R  is  gas  constant,  A  is  pre-exponential,  T  is 
temperature  and  fl  is  dimensionless  number  of  order  one. 

Hence,  Arrhenius  parameters  of  reaction  (9)  is  written  as: 


,  /q,  A,  exp  ( Jn) 

*■-£-  we  =A2exp 


—T„ 


(16) 


Viscosity,  thermal  conductivity  and  specific  heat  are  calculated 
by  considering  the  average  of  species  mass  fraction.  Radiation  heat 
transfer  between  micro-combustor  walls  is  considered.  Discrete 
ordinates  (DO)  model  is  used  to  simulate  surface-to-surface  radia¬ 
tion  heat  transfer.  Indeed,  non-slip  internal  walls  are  assumed  and 
no  species  flux  normal  to  the  surfaces  is  considered.  It  is  assumed 
that  heat  losses  from  micro-combustor  surfaces  to  the  ambient  by 
radiation  and  convection  heat  transfer. 


q  =  h(Ts<0-T00)  +  ec 7(^0 -O 


(17) 


Eddy  dissipation  concept  (EDC)  is  used  to  model  turbulence- 
chemistry  interaction.  EDC  includes  detail  of  chemical  mechanisms 
in  turbulent  flows  when  reaction  takes  place  in  small  turbulent 
conditions  or  so-called  fine  scales.  The  model  of  fine  scales  length 
fraction  is  modeled  br  Eq.  (18)  [48], 


(18) 


where  *  denotes  fine-scale  quantities. 

It  is  assumed  that  the  combustion  at  fine  scales  occurs  at  con¬ 
stant  pressure  when  the  current  species  and  temperature  in  the 
cell  is  taken  the  initial  conditions.  Indeed  Arrhenius  equation  gov¬ 
erns  the  speed  of  the  reactions  over  the  time  scale  presented  in  Eq. 
(19). 


(19) 


In  most  of  the  previous  investigations  about  macro-scale  flame¬ 
less  combustion  technology,  dilution  of  oxidizer  has  been  men¬ 
tioned  as  one  of  the  fundamentals  of  flameless  regime  formation. 
It  has  been  stipulated  that  when  the  percentage  of  oxygen  in  the 
combustion  air  increases  up  to  15%,  flame  is  constituted  [39],  The 
boundary  conditions  of  this  study  have  been  selected  based  on 
the  previous  experiments  in  macro-scale  flameless  mode.  In  con¬ 
ventional  combustion  the  temperature  of  oxidizer  (consist  of  21% 
02  and  79%  N2  by  vol.)  is  set  300  K  and  in  micro-flameless  mode 
the  temperature  of  inlet  oxidizer  (easel:  5%  02  and  95%  N2,  case2: 
7%  02  and  93%  N2  by  vol.)  is  adopted  900  K,  higher  than  self-igni¬ 
tion  temperature  of  methane. 

In  the  inlet  of  the  micro-flameless  combustor,  it  is  assumed  that 
u  =  cte,  T  =  cte  and  (T)  =  —  0. 

Indeed,  in  the  exhaust  %  =  ^  =  0  and  %  =  %  =  A 

*  ox  ay  az  ax  ay  az 

summary  of  the  boundary  conditions  of  the  simulation  is  pre¬ 
sented  in  Table  1. 


3.  Results  and  discussion 

The  grid  independent  test  is  done  by  changing  the  number  of 
elements  to  the  finer  elements.  The  number  of  meshes  were  set 
206,190,  254,202  and  312,021  respectively  to  observe  the  changes 
of  the  results.  However,  meaningful  variations  were  not  observed 
in  the  recorded  results  therefore  the  numerical  computation  is 
done  with  40,073  nodes  and  206,190  elements  to  reduce  the  num¬ 
ber  of  iterations.  The  validation  of  the  numerical  records  of  micro¬ 
combustion  in  conventional  mode  was  done  by  comparing  the 


124 


S.E.  Hosseini,  M.A.  Wahid/Energy  Conversion  and  Management  88  (2014)  120-128 


Table  1 

Boundry  condition  of  simulation. 


Properties 

Premixed  conventional  micro-combustion 

Micro-flameless  combustion 

Velocity 

Uniform  14*,  mixture  =  10  m/s 

Uniform  Viniet  oxidizer  =  30  m/s 

Violet  CH4=  1.32  m/s 

Temperature 

Uniform 

Uniform 

T'inlet  mixture  =  200  K 

^inlet  oxidizer  =  900  K 

Tinlet  CH4  =  300  K 

Oxidizer  ingredients 

N2:  93-95% 

02:  5-7% 

Density  of  oxidizer  (kg/m3) 

- 

0.382-0.383 

Density  of  mixture/fuel  (kg/m3) 

1.238 

0.770 

Pressure  gage 

OPa 

OPa 

Wall  material 

Steel 

Steel 

Wall  slip 

Non-slip 

Non-slip 

Thermal  conduction 

Mixed  (convection  +  radiation) 

Mixed  (convection  +  radiation) 

results  with  previous  experimental  [49]  and  numerical  [21  ]  results. 
Fig.  2  compares  the  axial  and  radial  temperature  profile  of  simu¬ 
lated  micro-flameless  mode  with  respect  to  various  numbers  of 
meshes. 

3.1.  Premixed  and  non-premixed  conventional  micro-combustion 

Most  of  previous  investigations  about  conventional  micro-com¬ 
bustion  were  done  in  premixed  circumstances  where  fuel  and  oxi¬ 
dizer  were  well  stirred  before  injection  to  the  micro-combustor. 
Due  to  the  small  length  and  diameter  of  micro-combustor,  applica¬ 
tion  of  premixed  combustion  method  could  be  the  best  strategy  to 
optimize  the  efficiency  of  the  combustor.  Moreover,  the  stability  of 
the  flame  and  efficiency  of  micro-combustor  could  be  enhanced  by 
using  bluff-body  [43],  In  this  regard,  a  triangular  bluff-body  is 
designed  at  the  inlet  of  the  micro-combustor  to  notice  the  effect  of 
bluff-body.  Depicts  temperature  distribution  inside  the  micro-com¬ 
bustor  in  premixed  and  non-premixed  conventional  combustion. 

Based  on  the  concept  of  chemical  kinetics,  chemical  reaction 
time  reduces  when  the  reaction  rate  increases.  In  the  micro-com¬ 
bustors,  higher  reaction  temperature  ensures  reduction  of  reaction 
time.  Hence,  enhancement  of  flame  temperature  as  well  as  heat 
loss  mitigation  from  the  micro-combustor  could  ensure  the  flame 
stability  50],  From  Fig.  3,  it  can  be  construed  that  making  stable 
flame  in  the  non-premixed  conventional  micro-combustion  is  not 
as  convenience  as  premixed  combustion  even  when  bluff-body  is 
employed.  In  the  other  hand,  in  premixed  conventional  micro¬ 
combustion  flame  is  constituted  in  the  downstream  region  of  the 
combustor  and  it  makes  some  limitation  for  inlet  velocity  of  the 
mixture,  thus  it  is  preferred  to  apply  premixed  combustion  with 
bluff-body  in  conventional  combustion  of  micro-scales. 

3.2.  Micro-flameless  combustion 

As  it  was  mentioned,  in  flameless  regime  high  temperature 
diluted  oxidizer  injected  to  the  furnace,  hence  for  flameless  com¬ 
bustion  mode  non-premixed  regime  is  utilized.  In  order  to  enhance 
the  stability  of  flameless  mode,  a  bluff-body  is  used  in  the  inlet  of 
the  micro  combustor.  Fig.  4  displays  the  structure  of  flameless 
micro-combustion  in  various  cases.  Indeed,  temperature  distribu¬ 
tion  along  the  centerline  and  wall  temperature  of  different  condi¬ 
tions  of  micro-flameless  mode  has  been  demonstrated. 

Micro-flameless  combustion  is  demonstrated  better  perfor¬ 
mance  in  terms  of  flameless  stability  and  temperature  distribution 
when  bluff-body  is  applied  at  the  inlet  of  combustor.  Fig.  4(A)  and 
(B)  illustrate  that  the  peak  of  temperature  shifted  to  the  inlet  of  the 
micro-flameless  combustor  (5%  and  1%  oxygen  concentration) 
when  the  bluff-body  was  employed.  The  maximum  temperature 
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Fig.  2.  (A)  Axial  and  (B)  radial  temperature  profile  of  simulation  (at  X  =  0.0075  mm 
from  the  inlet). 


in  bluff-body  micro-flameless  combustion  of  CH4  was  recorded 
1551  K  and  1376  K  when  oxygen  concentration  in  the  oxidizer 
was  7%  and  5%  by  volume  respectively. 

From  previous  studies,  it  was  found  that  both  C02  and  N2  are  eli¬ 
gible  to  be  applied  as  diluent  to  dilute  the  oxidizer  in  flameless 
mode  [51  ].  In  order  to  evaluate  the  performance  of  micro-flameless 
mode  in  different  conditions,  two  various  cases  were  solved  (case 
one:  02  concentration  is  7%  and  N2  is  applied  as  diluent,  case  two: 
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Premixed  conventional  micro-combustion  Non-premixed  Conventional  micro-combustion 


(A) 


Axial  distance  (mm) 


(B) 


Fig.  3.  (A)  Flame  structure  in  premixed  and  non-premixed  conventional  micro-combustion  with  and  without  bluff-body  and  (B)  temperature  distribution  along  centerline  of 
the  chamber. 


I 

i 


1.5130+003 
1.4376+003 
1.3816+003 
1.2860+003 
1.2106+003 
1.1346+003 
1  058O+003 
9.824O+002 
9.0660+002 
8.3070+002 

679 10+002 
6.0330+002 
5.275e+002 
4.5160+002 
3.7580+002 
3.0000+002 


[K] 


Micro-flameless  mode,  5%  Oxygen  concentration,  with  and  without  bluff-body 


Micro-flameless  mode,  7%  Oxygen  concentration,  with  and  without  bluff-body 
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Fig.  4.  (A)  Temperature  contour  of  micro-flameless  combustion,  (B)  temperature  distribution  along  the  centerline  and  (C)  wall  temperature  of  the  micro-flameless  mode. 
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Oxygen  concentration  is  7%,  diluted  by  N2 


Oxygen  concentration  is  7%,  diluted  by  C02 


Fig.  5.  The  effects  of  various  diluent  on  the  temperature  distribution  with  in  the 
micro-flameless  mode. 


— • —  Premixed  conventional  flame  with  bluff-body 
o  Flameless  combustion  with  bluff-body  5%  oxygen  concentration 
— Flameless  combustion  with  bluff-body  7%  oxygen  concentration 
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Fig.  6.  (A)  The  outlet  velocity  profile  and  (B)  the  outlet  temperature  profile. 


02  concentration  is  7%  and  C02  is  applied  as  diluent).  The  maximum 
temperature  of  combustor  was  recorded  1551  K  and  1345  K  when 
N2  and  C02  were  applied  as  diluent  respectively.  Therefore,  when 
N2  is  applied  as  diluent,  micro-flameless  system  shows  better  per¬ 
formance  in  terms  of  achieving  higher  temperatures.  The  lower 
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Fig.  7.  Wall  temperature. 


temperature  in  case  two  could  be  attributed  to  the  high  radiation 
heat  transfer  via  C02  particles  in  the  micro-combustor.  Fig.  5  illus¬ 
trates  the  effects  of  various  diluent  on  the  temperature  distribution 
with  in  the  micro-flameless  mode. 

In  order  to  interpret  the  micro-combustion  phenomena  in  both 
conventional  and  flameless  mode,  Damkohler  (D0)  number  should 
be  taken  into  account.  The  interaction  between  mixing  (or  flow 
time  scale)  and  reaction  in  a  combustion  system  is  represented 
by  Da  number  [39], 

Da  =  Xf/Xc  (20) 

where  z /  is  flow  time-scale  and  zc  is  chemical  time-scale. 

Large  Da  numbers  occurs  in  mixing  controlled  flames  and  it 
depicts  that  reactants  and  products  are  quickly  stirred  in  a  high 
turbulence  circumstance,  in  contrast,  low  Da  numbers  demon¬ 
strates  slow  chemical  reactions  [52],  Based  on  the  calculations, 
both  reaction  rates  (1?!  from  Eq.  (8)  and  R2  from  Eq.  (9))  in  the 
micro-flameless  combustion  mode  with  7%  02  concentration 
(Ri  =  0.7433,  R2  =  0.7195)  is  higher  than  conventional  micro-com¬ 
bustion  (R]  =  0.6915,  R2  =  0.6335).  In  addition,  Da  number  is  much 
greater  in  micro-flameless  regime  and  thus,  complete  and  stabile 
energy  conversion  is  expected  in  micro-flameless  combustion 
mode. 

In  conventional  premixed  micro-combustion,  when  the  mixture 
inlet  flow  rate  increases,  the  flame  length  becomes  longer  and  the 
flame  is  formed  further  downstream  of  the  inlet  of  the  combustor. 
Therefore,  longer  heating  length  is  needed  for  high  flow  rate  mix¬ 
tures  to  be  heated  up  and  reach  ignition  temperature.  Hence,  the 
velocity  of  the  mixture  is  one  of  the  limitations  of  conventional 
micro-combustors.  Since  ignition  is  eliminated  in  flameless  mode 
and  visible  flame  is  not  constituted,  higher  velocities  could  be 
adopted  in  micro-flameless  mode,  which  could  be  ideal  condition 
for  micro-turbines.  Furthermore,  the  temperature  of  flue  gases 
influences  significantly  on  the  efficiency  of  gas  turbine.  It  means 
that  high  exhaust  gas  temperatures  could  enhance  the  efficiency 
of  the  micro-turbine  [53],  Fig.  6  demonstrates  the  outlet  velocity 
and  temperature  profile  of  micro-flameless  regime  in  the  various 
cases. 

Based  on  Fig.  7,  the  results  show  that  the  outlet  velocity  as  well 
as  outlet  temperature  of  exhaust  gases  of  micro-flameless  combus¬ 
tion  (7%  oxygen  concentration)  is  significantly  higher  than  micro¬ 
flameless  mode  with  5%  oxygen  concentration  and  conventional 
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micro-combustion.  In  Fig.  3  it  was  shown  that  the  maximum  tem¬ 
perature  of  conventional  premixed  micro-combustion  is  around 
2200  K,  which  depicts  flame  temperature.  However,  from  Fig.  6B 
outlet  temperature  of  exhaust  gases  is  around  1300  K.  In  contrast, 
maximum  temperature  and  outlet  temperature  of  flue  gases  in 
flameless  mode  are  recorded  1520  K  and  1500  K  respectively.  It 
means,  temperature  inside  the  micro-combustor  is  more  homoge¬ 
neous  in  flameless  mode  and  hot  spots  are  eliminated.  Therefore, 
the  body  of  micro-combustor  could  be  protected  from  elevated 
temperatures.  Fig.  7  shows  the  wall  temperature  in  conventional 
and  flameless  mode.  The  fluctuation  of  wall  temperature  in  con¬ 
ventional  mode  ruins  the  material  of  the  combustor  wall.  Indeed, 
these  fluctuations  have  negative  effects  on  the  emitter  efficiency 
of  the  micro-combustor  in  micro-TPV  application. 

The  other  important  characteristic  of  micro-flameless  mode  is  its 
low  fuel  and  oxidizer  consumption.  Since  micro-combustors  are 
usually  applied  in  the  sensitive  situations,  it  is  vital  to  decrease  the 
size  of  fuel  and  oxidizer  cylinders  and  increase  the  period  of  their 
recharge  simultaneously.  The  rate  of  CH4-oxidizer  consumption  in 
premixed  conventional  micro-combustion  and  non-premixed 
micro-flameless  was  recorded  2.77  x  1CT3  g/s-2.63  x  1CT2  g/s  and 
8.05  x  10~4g/s-2.7  x  10~2  g/s  respectively.  It  means  that  with 
respect  to  the  specific  energy  generation,  consumption  of  a  specific 
volume  of  CH4-oxidizer  in  micro-flameless  takes  longer  time  and  the 
period  of  cylinder  recharging  is  prolonged. 


4.  Conclusion 

The  characteristics  of  micro-flameless  combustion  was  investi¬ 
gated  by  3D  simulation.  The  effect  of  bluff-body  on  the  tempera¬ 
ture  distribution  and  blow-off  limit  of  micro-flameless  mode  was 
simulated.  The  results  confirm  that  the  stability  of  the  premixed 
conventional  micro-flame  increases  when  a  triangular  bluff-body 
is  applied.  Indeed,  combustion  is  more  stable  in  micro-flameless 
mode  when  bluff-body  is  used.  Micro-flameless  combustion  with 
bluff-body  and  IX  02  concentration  showed  better  performance 
(in  terms  of  higher  peak  of  temperature  and  higher  flue  gases 
velocity)  than  other  cases  in  the  flameless  modes.  Moreover,  when 
N2  is  used  as  diluent,  the  maximum  temperature  of  micro-flame- 
less  mode  increases  compared  to  C02  utilization.  It  was  demon¬ 
strated  that  the  peak  temperature  of  premixed  conventional 
micro-combustion  (around  2200  K)  is  much  higher  than  the  max¬ 
imum  temperature  of  micro-flameless  combustion  (around 
1520  I<).  However,  the  temperature  of  exhaust  gases  (which  plays 
crucial  role  on  the  performance  of  the  micro-turbines)  in  conven¬ 
tional  mode  (1300  K)  is  lower  than  micro-flameless  regime 
(1500  I<).  The  fluctuations  of  temperature  in  the  conventional  com¬ 
bustor  wall  could  damage  the  inner  wall  of  the  micro-combustor, 
but  in  the  flameless  mode,  the  wall  temperature  is  moderate  and 
uniform.  Furthermore,  these  fluctuations  have  negative  effects  on 
the  emitter  efficiency  the  performance  of  micro-TPV  could  be 
reduced.  With  respect  to  the  specific  energy  generation,  the  rate 
of  CH4/oxidizer  consumption  in  micro-flameless  mode  is  lower 
than  conventional  combustion.  Some  advantages  of  micro-flame- 
less  combustion  could  be  summarized  as  below: 

•  In  the  micro-flameless  system,  combustion  stability  and  tem¬ 
perature  uniformity  inside  the  micro-combustor  increase  com¬ 
pared  to  the  conventional  mode. 

•  Reduction  of  fuel  and  oxidizer  consumption  could  be  obtained 
in  micro-flameless  mode,  thus  the  combustion  system  package 
is  more  compact. 

•  The  flame  quenching  which  has  been  mentioned  as  one  of  the 
most  important  limitations  of  the  conventional  micro-combus¬ 
tion  could  be  eliminated  in  the  micro-flameless  mode. 


•  Due  to  lift-off  of  the  flame,  the  inlet  velocity  of  the  mixture  is 
limited  in  the  conventional  micro-combustion,  which  could  be 
a  weak  point  in  the  micro-turbine  systems.  However,  in  the 
micro-flameless,  higher  inlet  velocities  could  be  achieved  due 
to  the  elimination  of  the  flame. 

•  In  micro-flameless  combustion,  inside  temperature  of  the  com¬ 
bustor  is  lower  than  conventional  mode.  Hence,  a  broader  range 
of  materials  could  be  employed  in  a  micro-flameless  system. 
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